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ABSTRACT 

Gamma-ray bursts (GRBs) are tremendous explosions visible across most of 
the Universe, certainly out to redshifts of z = 4.5 and likely out to 2; ~ 10. 
Recently, GRBs have been found to have a roughly constant explosive energy as 
well as to have two luminosity indicators (the spectral lag time and the variability) 
that can be used to derive the burst's luminosity distance from the gamma- 
ray light curve alone. There currently exists enough information to calibrate 
luminosity distances and independent redshifts for nine bursts. From these, a 
GRB Hubble diagram can be constructed, where the observed shape of the curve 
provides a record of the expansion history of our Universe. The current nine 
burst diagram is sparse, yet formal limits can be placed on the mass density 
of a flat Universe. This first GRB Hubble diagram provides a proof of concept 
for a new technique in cosmology at very high redshifts. With the launch of 
the SWIFT satelhte in 2003, we should get ~ 120 bursts to produce a Hubble 
diagram impervious to all effects of dust extinction and out to redshifts impossible 
to reach by any other method. 

Subject headings: cosmology: cosmological parameters — cosmology: distance 
scale — cosmology: observations — gamma rays: bursts 



1. Introduction 

As always in astronomy, the determination of distances is a crucial and difficult problem. 
Until recently, the distance scale of GRBs was unknown by over 12 orders-of-magnitude. In 
1997, the discovery of optical and radio counterparts (Costa et al. 1997; van Paradijs et al. 
1997; Frail et al. 1997) proved that at least the long-duration bursters were at cosmological 
distances with redshifts of z ~ 1. The measurement of GRB redshifts requires deep optical 
spectra, and to date only 24 redshifts are known for bursts with unknown selection effects. 

If GRB distance indicators can be found that use only gamma-ray data, then we can 
measure the demographics and cosmology of large and well-understood samples of bursts. 
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Two such luminosity (and hence distance) indicators have recently been proposed. The 
first (Norris, Marani, & Bonnell 2000) relates the burst luminosity (L) with the spectral 
lag (riag), which can be idealized as the time between peaks as recorded at high and low 
photon energies. The second (Fenimore & Ramirez-Ruiz 2000) relates the luminosity with 
the variability {V), which is a specific measure of the 'spikiness' of the burst light curve. 
High luminosity bursts have short lags and spiky light curves, while low luminosity events 
have long lags and smooth hght curves. These two relations make GRBs into 'standard 
candles', in the same sense as for Cepheids and supernovae where an observed light curve 
property can yield the luminosity and then the distance. 

The two luminosity indicators were originally proposed and calibrated with 6 or 7 bursts. 
The addition of 3 or 2 further bursts (for a total of nine) have fallen on the original relations, 
hence adding confidence in their utility. Further, if both the L/riag and L/V relations are 
true, then there must be a particular Tiag/V relation, and this prediction has been strongly 
confirmed with an independent sample of 112 BATSE bursts (Schacfcr, Deng, & Band 2001). 
Finally, the very long lag bursts have been shown to be very low luminosity as demonstrated 
by their concentration to the local supergalactic plane (Norris 2002). 

These luminosity indicators have been used to identify specific bursts (Fenimore & 
Ramirez- Ruiz 2000) that arc at redshifts of z ~ 10 as well as to show that the star formation 
rate of the Universe is rising steadily (Fenimore & Ramirez-Ruiz 2000; Schacfer, Deng, & 
Band 2001; Lloyd- Ronning, Fryer, & Ramirez-Ruiz 2001) from 2; ~ 2 to 2; > 6. This Letter 
reports on the construction of a Hubble diagram (a plot of luminosity distance, D^, versus 
redshift) as a means of measuring the expansion history of our Universe. 

2. GRB Hubble Diagram 

Only nine GRBs have the required information of red shift (2;), peak flux (P), lag time 
(Tiag)i and variability {V). These data are collected in Table 1, along with the characteristic 
photon energy (Epeak) and the observed luminosity {Lobs)- These nine bursts were all de- 
tected by BATSE with redshifts measured from optical spectra of either the afterglow or the 
host galaxy. The highly unusual GRB980425 (associated with supernova SN1998bw) is not 
included because it is likely to be quahtatively different from the classical GRBs. Bursts with 
red shifts that were not recorded by BATSE cannot (yet) have their observed parameters 
converted to energies and fluxes that are comparable with BATSE data. 

Simplistically, plots of Lobs versus Tiag and Lobs versus V can calibrate the luminosity 
indicators, which then can yield luminosity distances to each burst for plotting on a Hubble 
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Diagram and fitting to cosmological models. In practice, there must be a simultaneous chi- 
square minimization for the luminosity calibrations and the cosmology, so as to avoid the 
effects of possible correlations between the indicators and distance. 

If the Earth is not along the central axis of the GRB's jet, then various off- axis effects 
will change the observed properties. For light viewed off- axis by angle 6 from a jet moving 
at j3 times the speed of light, the transverse Dopplcr shift, redshift, and beaming will change 
the various observed properties by powers oi B = {1 — (3) / {1 — j3cos9). The observed value 
of Epf^ak will scale as -B, Tiag will scale as -B"^, V will scale as the inverse of a time and 
hence as i?, while L will have a model dependent variation that scales roughly as . For 
structured jets (Rossi, Lazzati, & Rees 2002; Zhang & Mcszaros 2002), the luminosity and 
bulk Lorentz factor can be parameterized as being proportional to some positive power of 
B. When viewed on-axis, the bursts will presumably display fairly tight lag/luminosity 
and variability/luminosity relations. These relations can be converted to similar relations 
involving observed off-axis quantities with an extra factor of B to some power. Fortunately, 
the observed distribution of Ep^ak is remarkably narrow for bursts of a given brightness 
(Mallozzi et al. 1995; Brainard et al. 1998), and the removal of kinematic effects shows that 
the on-axis -Epeafc is virtually constant (Schaefer 2002). Thus, Ep^ak should be proportional 
only to some power of B, and this relation can be used to convert the B dependency of 
the luminosity relations into a Ep^ak dependency. Essentially, we can use the observed Ep^ak 
value to get information about the off-axis angle and correct the observed Tiag and V values 
to those which would be seen on-axis. In practice, the correction factor can be taken as 
some power of £'peafe(l + z)/AQ{)keV. (The division by 400 keV is to minimize correlations 
between the normalization constant and the exponent during the fits.) The exponents of the 
correction factors will be free parameters which depend on the scenario and jet structure. 

With the off- axis correction, the luminosities can be derived by fitting the relations 

L 0^ ^laZorr ^^d L OC V^^, whcrC Tiag,corr = Tiag X [Ep^aki'^ + z)/400]^'-« and Korr = X 

[Epeaki)- + -2)/400]'^^. Plots of Lobs versus Tiag,corr and Lobs versus Vcorr for the nine bursts 
with known redshifts (see Figure 1) have slopes of aiag and ay, while the scatter in these 
plots will be minimized for the best values of eiag and ey. This gives aiag = —1.27 ± 0.20, 
av = 1.57 ± 0.17, eiag = 0.6 ± 0.7, and ey = 0.85 ± 0.40. Thus: 

L = 10™-°3[(t,„, X {Ep,,,{l + z)/mr-T'-'': (1) 

L = 10''-''[{V X {E,eafe(l + ^)/400}°-«Y-'', (2) 

where L is in erg s~^, Tiag is in seconds, V is dimensionless, and Ep^ak is in keV. As in- 
dependent checks, a plot of Tiag^corr versus Vcorr for 93 BATSE bursts will have a slope of 
Oiiag/otv-i and the scatter will be minimized for the best values of eiag and ey. This gives 
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OLiagjoLy — —0.7 ± 0.2, eiag — 0.0 ± 0.3, and ey — 1.1 ± 0.3; values which are consistent with 
the cahbration from Figure 1. 

The lag and the variability produce two independent luminosity values through Eqs 1 
and 2. These can be combined as weighted averages (with the weights equal to the inverse 
square of the one-sigma errors) to produce a combined luminosity {Lcomh) and then Dl, as 
given in Table 1. A plot oi Dl versus z is a Hubble diagram (Figure 2). 

Equations 1 and 2 are optimal for studies of GRBs where the cosmology is not in ques- 
tion. When the goal is to use the bursts as cosmological markers, care must be taken to avoid 
circular logic. In particular, the previous paragraph assumed a specific cosmology (Q^ = 0.3 
for a flat Universe and Hq = 65 km Mpc~^) for calculating the Dl values as a function 
of z. To avoid this circularity, the fits must be performed with the cosmological parameters 
as additional fit parameters. That is, the model being fit to the data will have cosmological 
plus GRB parameters that are to be simultaneously determined by the minimization of the 
chi-square. This will tie the cosmology conclusions into GRB calibrations, but both will be 
well constrained when many bursts are available. A correct procedure to find the best value 
of in a flat Universe is (1) to flx ft^, (2) derive and L^hs for each burst for that 
cosmology, (3) flt Lg^g as power laws of Tiag^corr and Vcorr, (4) use the best flts to derive Lcomb 
and then Dl for each burst, (5) calculate chi-squarc by comparing the model and observed 
Dl values, (6) repeat steps 1-5 to identify the best fit Qm value as that where chi-square 
is minimized, and (7) identify the onc-sigma range where chi-squarc is within unity of its 
minimum value. By this procedure, the lowest chi-squarc corresponds to Qm = 0.05 with 
the one-sigma constraint < 0.35. The three-sigma limit is larger than unity, so this 
particular measure is not helpful for cosmology. This new result is completely independent 
of, yet in agreement with, those from supernovae and other methods (Durrer & Novosyadlyj 
2001). However, the current value for Q^n is not yet of high rehabihty because the number 
of degree of freedom in the fit are small and the slope of the Hubble diagram is dominated 
by just two high redshift events. While this new result does not seriously constrain current 
cosmology, it does act as a demonstration of principle and a sign for the future. 



3. Comparison with Supernovae 

A comparison with the supernova-based Hubble diagram is inevitable. GRBs have both 
advantages and disadvantages when compared to supernovae. 

At least currently, supernovae have an advantage over GRBs for the accuracy of the 
distance measurements for a single event. Well-observed nearby Type la supernovae have 
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an RMS scatter about their Hubble Diagram of 0.18 mag (Phillips et al. 1999) which 

translates into an uncertainty in the log of the distance of 0.04. GRBs currently have an 
RMS scatter about their calibration curves of 0.35 and 0.20 in the logs of lag and variability; 
which translates into uncertainties in the log of the distance of 0.22 and 0.16 respectively. 
These two independent measures can be combined to give an uncertainty of 0.13 in the log of 
the distance. Thus, individual supernovae are roughly three times more accurate as standard 
candles than are individual GRBs. 

The physics of gamma-ray emission from relativistic shocks in GRBs is largely known 
(but not how to create the magnetic fields) although the explosion scenario is still uncer- 
tain. The basic physics of supernovae is well known (although the Type la progenitors are 
still not identified). Supernovae can have their decline-rate versus luminosity calibration 
from low-redshift events, and this allows the elegance of being able to make the calibration 
substantially separate from the cosmology. Thus, supernovae are much better understood 
than are GRBs, and this will be a substantial comfort for using them as standard candles. 
Nevertheless, it is unclear where this advantage for supernovae will pay off. Currently, the 
analyses of both supernova and GRB Hubble Diagrams are entirely empirical with no con- 
tribution from theory, so a better theoretical understanding of supernovae has not helped 
in any specific manner. Supernova will have better known evolution effects; yet Branch et 
al. (2001) show that evolution is no problem for supernovae and their argument applies 
identically to GRBs. 

Supernovae have substantial problems with a higher dust density at high red shift 
(Totani & Kobayashi 1999) and the possibility of grey dust (Aguirre 1999). Gamma 
radiation suffers from no extinction. 

After a tremendous observational effort, the current supernova Hubble diagram (Perl- 
mutter et al. 1999) extends only out to z = 0.97. (SN1997ff at z 1.7 has very large 
uncertainties (Riess et al. 2001) even without the large corrections for gravitational lensing 
(Lewis & Ibata 2001; Moertsell, Gunnarsson, & Goobar 2001) and the uncertainty of the 
supernova type.) The dedicated SNAP satellite (SNAP collaboration 2002), proposed for 
launch in 2008, will go out to 2; = 1.7 with exquisite accuracy in the light curves. In contrast, 
the GRB Hubble Diagram is already out to z = 4.5 and likely it will be extended to 2; ~ 10 
(Fenimore & Ramirez-Ruiz 2000; Lloyd-Ronning, Fryer, & Ramirez-Ruiz 2001) or farther 
(Lamb & Reichart 2000; Bromm & Loeb 2002). With the launch of the SWIFT satellite 
in 2003, the GRB Hubble Diagram will be available for the redshift ranges 1 < z < 1.7 and 
l-"^ < ^ ^ 10 roughly five years before SNAP can extend the supernova Hubble Diagram 
only to 1 < z < 1.7. 

In summary, GRBs are ~ 3x worse in accuracy than supernovae, but this is traded- 
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off for the lack of problems from extinction, extension to z ~ 10, and an answer many 
years before SNAP. I do not think that the GRB Hubble Diagram will replace the supernova 
Hubble Diagram. The reason is that hoth diagrams will have known and unknown systematic 
problems which will make results from any one method not conclusive. What is needed is 
the concurrence of multiple independent methods. So hoth GRB and supernova Hubble 
Diagrams are needed for a confident result. 

4. Implications 

Conceptually, the biggest question related to a construction of the GRB Hubble dia- 
gram is the possibility of luminosity evolution. Fortunately, there are two strong reasons that 
any such effects must be small. First, theoretical explanations for the luminosity indicators 
(Schaefer 2001; loka & Nakamura 2001; Plaga 2001) tie the relations to simple physics (in- 
volving relativistic effects operating on any light source, the delay in light travelling different 
paths, and the cooling of any body by radiation) that are unrelated to properties (such as 
metallicity) that might evolve. Second, just as for the purely empirical supernova Hubble 
diagram (Branch et al. 2001), any drift in population average properties (Lloyd- Ronning, 
Fryer, & Ramirez-Ruiz 2001) is irrelevant as the calibrations will work for any individual 
event whether near or far. That is, it does not matter whether population average properties 
(like luminosity) evolve because the distances are derived for individual events that are each 
correctly calibrated. 

Within two years, the GRB situation will be greatly improved with the dedicated SWIFT 
satellite (Gehrels 2000) currently scheduled for launch in September 2003. SWIFT is ex- 
pected to get ~ 100 bursts per year, of which ~ 40% will have directly measured redshifts. 
The redshift will be difficult to measure for 2; > 6 bursts, yet it is possible by locating the 
Lyman break with broad-band photometry or spectroscopy extending into the near infrared. 
In its nominal three-year lifetime, we can get ~ 120 bursts with measured values for 2;, P, 
'^lagi V, a-nd Epgak- This many bursts will allow for high accuracy calibration of the lumi- 
nosity indicators. The large number of bursts will also allow for the reduction of statistical 
errors (dominated by the intrinsic scatter of the bursts) by a factor of ~ 10. If the bursts are 
divided into six redshift bins, each bin will have an error of 0.03 in the log of the distance. 
This accuracy will, for example, yield an uncertainty of ~ 0.03 in the derived value of Qm for 
a flat Universe, even with no improvements in the luminosity indicators. So with SWIFT, 
the prospect is that we can produce an accurate GRB Hubble diagram with ~ 120 bursts 
from O.K ^ < 10 by 2006. 

What can we learn from the GRB Hubble Diagram? We can test the predicted shift 
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of the Universe from matter to dark-energy domination in the range 1 < z < 2. We will 
also look for various predicted quintessence (Weller & Albrecht 1999) and non-standard 
(Mannheim 2002) effects in the range 2 < z < 10. The typical sizes of these effects are 0.4 
in the log of the distance around a redshift of z = 3. The current paradigm of cosmology 
(the new inflation perhaps with quintessence) is so new and untested that surprises could 
easily await in the unknown regime of 2; > 1. 
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Table 1. GRB Luminosity Distances From Lags and Variabilities. 



GRB 


z 


pa 


'lag 


yc 






Log{LcombY 


Dz,8 


970508 


0.84 


1.2 ±0.1 


0.307 ± 0.065 


0.0010 ±0.0010 


137 ±14 


50.89 ± 0.04 


50.70 ± 0.39 


4600 ± 2100 


970828 


0.96 


4.9 ±0.1 


0.028 ±0.007 


0.0078 ±0.0006 


176 ±4 


51.65 ±0.01 


51.97 ±0.26 


9600 ± 2900 


971214 


3.41 


2.3 ±0.11 


0.010 ±0.004 


0.0175 ±0.0012 


107 ±6 


52.69 ±0.02 


52.62 ±0.27 


29500 ± 9100 


980703 


0.97 


2.6 ±0.12 


0.147 ±0.056 


0.0025 ±0.0005 


181 ± 9 


51.38 ±0.02 


51.16 ±0.28 


5200 ± 1700 


990123 


1.60 


16.6 ± 0.24 


0.015 ±0.005 


0.0120 ± 0.0005 


267 ±3 


52.74 ±0.01 


52.49 ± 0.26 


9500 ± 2900 


990506 


1.30 


22.2 ±0.27 


0.011 ±0.004 


0.0320 ± 0.0080 


280 ± 30 


52.64 ±0.01 


52.94 ±0.29 


13800 ± 4600 


990510 


1.62 


10.2 ±0.2 


0.012 ±0.003 


0.0352 ±0.0014 


74 ±2 


52.54 ± 0.01 


52.65 ±0.26 


14600 ± 4400 


991216 


1.02 


82.1 ±0.5 


0.0050 ± 0.0020 


0.0152 ± 0.0003 


250 ±3 


52.95 ±0.01 


52.68 ±0.27 


5400 ± 1600 


000131 


4.5 


1.8 ±0.2 


0.0009 ± 0.0004 


0.0121 ±0.0018 


183 ± 15 


52.86 ± 0.05 


53.08 ±0.28 


57000 ± 18500 



^Peak flux for the brightest 256 ms time interval in photon s ^cm ^ for 50-300 keV (Paciesas et al. 1999). 

''Lag in seconds as calculated from the peak of the cross correlation between BATSE channels 3 and 1 using data brighter than 
10% of the peak flux (Norris, Marani, & BonncU 2000; Norris 2002). These lags have been corrected by 1 + z to the frame of 
the GRB. The quoted one-sigma measurement uncertainties in the log of the lag should be added in quadrature with 0.35, the 
population dispersion, to get the effective uncertainty for luminosity determination. 

''Variability following the definition of Fenimore & Ramirez- Ruiz (2000); where V is the normalized variance of the light curve 
around a O.lSTgo box-smoothed light curve. The quoted one-sigma measurement uncertainties in the log of the variability should 
be added in quadrature with 0.20, the population dispersion, to get the effective uncertainty for luminosity determination. 

''Photon energy (in kcV) of the peak in the i/F,^ spectrum as measured by Mallozzi ct al. (1995) with BATSE data. 

"The base-10 logarithm of the observed luminosity in erg s~^; calculated as 4nPeDj^, where e ~ 1.7 X lO^^erg as the average 
energy of a photon in the 50 — SOOkeV range and is the luminosity distance which is taken from z for an assumed flat Universe 
with Qrn = 0.3 (with equations 14, 15, 16, and 21 of Hogg 1999). The calculated observed luminosity depends on the adopted 
Hubble constant (here taJjen as 65 km s~^Mpc~^); although the cosmology only depends on the shape of the Hubble diagram and 
hence is independent of the Hubble constant. 

'The base-10 logarithm for the combined luminosity in erg s~'; calculated as the weighted average of the results of Eqs 1 and 
2 (for which Qm = 0.3 has been adopted). 

SThe luminosity distance in Mpc calculated as (Lcomb/^TrPe)"'®, for luminosity based on an Clm = 0.3 calibration. 
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Fig. 1. — Calibration curves for lag and variability. Plots of the observed luminosity versus 
lag and variability can be used to calibrate the luminosity indicators. Sufficient data are 
available for only 9 bursts with redshifts. The plots here can be fitted to produce Eqs 1 
and 2 (displayed as the straight lines in the two panels). The displayed error bars include 
not only the observed measurement uncertainties but also the intrinsic scatter for the burst 
population of 0.35 and 0.20 (for lag and variability respectively) in logarithmic units, as 
determined by making the reduced chi-square of the fits close to unity. 
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Fig. 2. — GRB Hubble Diagram. This Hubble diagram is for 9 bursts from 0.84 < z < 4.5. 
With SWIFT, we should soon have ~ 120 bursts from roughly 0.1 < 2; < 10, and this 
will allow a precision measure of the expansion history of the Universe to unprecedented 
high redshifts. The particular diagram above was made for luminosity indicators calibrated 
assuming a flat Vtm ~ 0.3 Universe (as in Figure 1 and Eqs 1 and 2) and the chi-square for the 
data versus the cosmological model (the curve in the figure) is 6.2. The best fit cosmology 
can be found by minimizing the resulting chi-square as the cosmological model varies both 
the calibration curves and the resulting Hubble diagram. 



